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We study neutrino energy emission rates (emissivities) due to electron 
bremsstrahlung produced by ee and ep collisions in the superfluid neutron 
star cores. The neutrino emission due to ee collisions is shown to be the 
dominant neutrino reaction at not too high temperatures (T < 10^ K) in 
dense matter if all other neutrino reactions involving nucleons are strongly 
suppressed by neutron and proton superfluidity. Simple practical expres- 
sions for the ee and ep neutrino emissivities are obtained. The efficiency of 
various neutrino reactions in the superfluid neutron-star cores is discussed 
for the cases of standard neutrino energy losses and the losses enhanced by 
the direct Urea process. 
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1. Introduction 

In the first 10^ years of their life, neutron stars (NSs) cool mainly via 
neutrino emission from their cores. This emission results from weak inter- 
action processes involving baryons and leptons. It is important, thus, to 
study various neutrino reactions. For simplicity, we will mainly consider 
the NS cores composed of neutrons n, with some admixture of protons p 
and electrons e. The neutrino reactions in the neutron-star cores are tradi- 
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tionally divided into two groups, which differ drasticahy by their efficiency 
and produce the standard or rapid cooUng (e.g., ref. [1]). 

The standard coohng is mainly provided by the neutron and proton 
branches of the modified Urea reactions [2, 3] 

n + N^p + N + e + Ue, p + e + N^n + N + Ue, (1) 

where is a nucleon (n or p), and by the nucleon-nucleon (nn, np, and pp) 
bremsstrahlung processes 

N + N N + N + u + . (2) 

The rapid cooling is strongly enhanced by the direct Urea processes 

n ^ p + e + Ue , p + e ^ n + Ue , (3) 

which can operate [4] only in the central regions of rather massive NSs with 
proton fraction exceeding 11% of the total baryon density. If the direct 
Urea reaction is allowed (in a non-super fluid NS), the neutrino emissivity is 
typically 4-5 orders of magnitude higher than in the standard reactions (|l]) 
and d). 

It is widely accepted (see, e.g., [5, 6] and references therein), that neu- 
trons and protons in a NS core can be in a superfluid state. The superfluidity 
is generally thought to be of BCS type, produced by attractive component 
of nuclear forces. Various microscopic theories predict very different crit- 
ical temperatures of neutron and proton superfluids. Ten and Tcp, which 
range from 10^ to 10^'^ K. The neutron superfluidity in the NS cores is 
generally believed to be produced by triplet-state pairing of neutrons (with 
anisotropic superfluid energy gap), while the proton superfluidity is due to 
the singlet-state pairing (with isotropic gap). If T ^ Ten and T <C Tcp the 
superfluidity strongly suppresses all the neutrino reactions (||)-(§) involving 
nucleons (e.g., [1, 3, 7]). 

Onset of the neutron or proton superfluidity with the decrease of temper- 
ature in a cooling NS initiates also an additional specific neutrino reaction 
which can be considered as neutrino production due to Cooper pairing of 
quasinucleons N _ _ 

N + N ^ u + i? . (4) 

The emission due to singlet-state pairing of neutrons has been calculated by 
Flowers et al. [8] and Voskresensky & Senatorov [9, 10] while the emission 
due to triplet-state pairing of neutrons and singlet-state pairing of protons 
has been determined by Yakovlev et al. [11, 12] and Kaminker et al. [13]. The 
effect is much stronger for neutrons than for protons due to the smallness 
of axial-vector electroweak currents for protons. 
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In contrast to the non-superfluid cores, where the main neutrino emis- 
sion is produced either by the modified or by the direct Urea processes 
(depending on the equation of state and density), very different neutrino 
mechanisms can dominate in the superfiuid cores, depending on the relative 
magnitudes of T, Ten and Tcp. In particular. Cooper-pair neutrino emission 
(§) can be very important if 0.1 Tcn < T^n (e.g., ref. [12]). At lower 
temperatures, T <C T^n^ this emission falls exponentially. 

On the other hand, as mentioned by Schaab et al. [14], an additional 
neutrino emission mechanism, the neutrino-pair bremsstrahlung due to ep 
scattering 

e + p^e+p + u + u, (5) 

can be important in the superfiuid NS cores. This mechanism, and the 
closely related mechanism of neutrino bremsstrahlung due to ee scattering 

e + e^e + e + iy + i> (6) 

has been considered briefiy by Kaminker et al. [15]. The latter authors have 
estimated the neutrino energy emission rates (emissivities) in these reac- 
tions using simplified criteria of similarity of neutrino emission and electron 
conduction processes. It is quite evident that these reactions, especially the 
neutrino ee bremsstrahlung, can be important in highly superfiuid NS cores, 
where all the reactions involving nucleons are drastically suppressed by the 
nucleon superfiuidity. Thus, rigorous calculation of the emissivities in the 
reactions (|5|) and (^) deserves special study. 

Notice that the neutrino bremsstrahlung due to ee scattering of nonde- 
generate electrons at not too high densities, p < 10^^ g cm~'^ (corresponding 
to the outer NS crust), was considered by Cazzola & Saggion [16, 17]. The 
authors obtained general but cumbersome expressions for the emissivity [16] 
and performed numerical calculations using Monte Carlo method [17]. The 
main result of their consideration was that the emissivity due to the ee 
bremsstrahlung did not exceed 1% of the total emissivity from other neu- 
trino reactions in a NS crust (pair annihilation, photoneutrino and plasmon 
decay). We will show that this is not so in a NS core. 

The aim of the present article is to carry out an accurate calculation of 
the neutrino emissivities for the bremsstrahlung reactions (^) and (|), and 
to analyze the efficiency of various neutrino production mechanisms (1 
in superfiuid NS cores. 
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2. Superfluid neutron-star cores 

Mass density of the matter in a NS core is expected to range from 
about 0.5po at its outer edge [18], to (5-10)/9o at the stehar center, where 
po = 2.8 X 10^^ g cm~^ is the saturation density of nuclear matter. At 
P ~ 2yOo, neutron-star matter consists mostly of neutrons, with a few per- 
cent admixture of protons, electrons, and possibly muons. The fraction of 
muons is usually much smaller than that of electrons, and a simple npe 
model is a good approximation. At higher densities, p > 2pQ, other parti- 
cles may appear, such as hyperons, condensed pions or kaons, or even free 
(deconfined) quarks. Neutrino emissivity of dense neutron-star matter with 
different compositions has been reviewed by Pethick [1]. For simplicity, we 
will neglect possible presence of hyperons and exotic particles, and consider 
npe matter only. All constituents of the matter are strongly degenerate. The 
electrons form an ultrarelativistic and almost ideal gas, with the electron 
Fermi- momentum 

PFe = HSTT^Ue)^/^ = 648.6 {ne/m)^^^ nieC , (7) 

where and me are electron number density and mass, respectively, and 
no = 0.16 fm~^ is the standard nuclear saturation density. The electron 
Fermi energy (chemical potential) is pe ~ PFeC ~ (100-300) MeV. Neutrons 
and protons are nonrelativistic and form strongly interacting Fermi liquids. 
The condition of beta-equilibrium implies that the chemical potentials of 
particles satisfy the equality /i„ = Pp + He- Charge neutrality leads to 
the equality of proton and electron number densities (up = n^) and Fermi 
momenta [p-pe = PFp)- The Fermi energy of neutrons is close to that of 
electrons, while the Fermi energy of protons is much smaller (typically, 3- 
30 MeV). 

Both nucleon species, n and p, can be in a superfluid state (Sect. 1). 
Microscopically, the superfluidity leads to the appearance of the energy gap 
A in the nucleon dispersion relation. For a triplet-state neutron superfluid, 
the gap is anisotropic along the neutron Fermi surface. For a singlet-state 
pairing of protons the gap is isotropic. The proton gap Ap = Ap(T) de- 
pends only on temperature T. In the frame of the BCS theory (which we 
adopt further on) , the proton critical temperature Tcp is related to the zero- 
temperature energy gap Ap(0) by Tcp = Ap(0)/(1.764 A;b) [19], where /cb is 
the Boltzmann constant. 

Superfluidity of neutrons and protons suppresses the main neutrino gen- 
eration mechanisms (|l|)-(^ in the NS core and opens an additional neutrino 
reaction (^). Strong superfluidity of neutrons and protons {Ten 3> T and 
Tcp ^ T) suppresses drastically all neutrino reactions (|l])-(|5|) involving nu- 
cleons. Superfluidity of protons affects even the neutrino ee bremsstrahlung 
(^), (not very strongly), by changing screening of ee interaction. 
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The plasma screening momentum Qs in npe matter is given by (e.g., ref. 
[20]) 

=2 



y. 



-i l + (B) 



Api Txfic \ m*p-pe 



where m* = /ie/c^ and m* is an effective proton mass (which determines 
the proton density of states near the Fermi level and which can differ from 
the bare mass nip due to polarization effects in dense matter). The first 
term in brackets comes from the electron screening, while the second term 
is due to the proton screening. The latter term contains the function Dp. 
For nonsuperfluid protons (T > T^p) one has Dp = 1 , while for T < T^p 
the factor Dp describes superfluid reduction of proton screening. Gnedin &: 
Yakovlev [20] fitted numerical values of this factor, calculated for singlet- 
state proton pairing, by 

1/2 



Dp = 0.9443 + a/(0.0557)2 + (0.1886y)^ 



X exp (^1.753 - ^(1.753)2 + y2j , (9) 

where y = lS.p{T) / [k-^T) . According to Levenfish & Yakovlev [21] the pa- 
rameter y at any T < T^p is given by the fitting expression 



|- ( 1.456 -0.157 ^1^ + 1.764 . (10) 



3. General formalism of neutrino ep and ee bremsstrahlung 

Consider two neutrino bremsstrahlung mechanisms (^) and (P) for stron- 
gly degenerate, relativistic electrons and strongly degenerate, nonrelativistic 
protons. We will treat these processes using the standard perturbation 
theory of weak interactions. We will use the units in which c = Ti = k-Q = 1 
and return to ordinary physical units in the final expressions. 

The emissivities Q^p and Qee [erg cm~^ s~^] of the processes (|5|) and (|6|) 
can be written as 

Q12 = j dPi j dp2 ydp'i j dp2 j dki. j dk'^ 

X <5(4)(p,+P2-P{-P^-K)a;/i/2(l-/()(l-/^)^, (11) 

where Pi = (ei,pi) and P2 = (e25P2) are 4-momenta of two charged parti- 
cles before scattering, P[ = (e'^, p'^) and P2 = (£2' P2) are 4-momenta of the 
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same particles after scattering. Index 1 labels an electron and index 2 labels 
either a proton or a second electron; (ep = 1, Cee = 1/2 (to avoid double 
counting of the same ee scattering events). The quantities K^^ = (cjj^jkjy) 
and Kl = (u;^,k'j,) are 4-momenta of neutrino and antineutrino, respec- 
tively, and K = Ky + = (w, k) is 4-momentum of the neutrino pair 
(ijj = Ui, + ul and k = kj^ + k^). Furthermore, 



fi 



1 + exp 



(12) 



is a Fermi-Dirac distribution of a particle before scattering {i=l, 2), and 
f- is the corresponding distribution for a particle after scattering. The 5- 
function describes 4-momentum conservation. Finally, W is the differential 
transition rate 



w = ^ 

2 (2u;,)(2c^^,)(2ei)(2e2)(2e'i)(2e'2) 



(13) 



where Gp = 1.436 x 10~^^ erg cm^ is the Fermi weak interaction constant, 
and |Mi2p is the squared transition amplitude. Summation is over spin 
states of charged particles before and after scattering a = (ai, a2, cr[,a'2) 
and over the neutrino flavors v = (i^e, f^, i^r)- The neutrino energies are 
assumed to be much lower than the intermediate boson mass (~80 GeV). 
The processes in question are represented by four {ep scattering ) or eight 
(ee scattering) Feynman diagrams. Each diagram includes a neutrino-pair 
emission four-tail vertex. The amplitudes are 



Mep = MA, Mee = MA-MB 



(14) 



where 



Ma = ^ne^l" 
1 



+ 



Mf 



hi - q\ 
47re2r 



{^'.I'u,) (.;l«u.) 

^ «l'u,) (4iS«0 



+ 



ti2 - qt 

L% = Tea Ge{P[ + i^) 7/3 + 7/3 Ge{Pl - K) 
41 = G2{P!, + iC) 7/3 + 7/3 G2{P2 - K)T 



(15) 



(16) 



2a 
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= Tea Ge{P2 + i^) 7/3 + 7/3 Ge{Pl " K) T^a , 
41 = Tea Ge{P{ + i^) 7/3 + 7/3 Ge(P2 " K) T^a , (17) 



Cey7" + CeA7V, = Cpv7" + CpA7 V , (18) 



= (Vii,7"(l + 'y^)'tp,y) is neutrino 4-current; ui, u'j^, U2 and U2 denote 
standard bispinor amplitudes of charged particles; ipi, and ipl are the stan- 
dard bispinor amplitudes of emitted neutrinos; upper bar denotes a Dirac 
conjugate, Greek indices a and /? run over (0,1,2,3), 7" and 7^ are Dirac ma- 
trices; Ge{P) and Gp{P) are the electron and proton propagators Gi{P) = 
{PqI" + ^i)/{P'^ — iTT-j)- We use the bare proton mass m2 = nip in the 
proton propagator and in the normahzation of the proton bispinor ampli- 
tudes UpUp = 2mp in Eqs. (^) and (p^). It will be shown that the final 
results are independent of nip which is a consequence of nonrelativistic ap- 
proximation for protons. Nevertheless the results do depend on the proton 
density of states near the Fermi level, i.e., on m*. An additional term Mb in 
Mee appears due to identity of colliding electrons; it contains four diagrams 
with interchanged final electron states. In Eqs. (|l^) and (16) we have used 



photon propagators in the Feynman gauge [22] and four familiar kinematic 
invariants 

tn = {Pi-P[f, t22 = {P2-P^f, 

t21 = {P2 - P[f, h2 = (Pi - P2?, (19) 

where tu and ti2 are the standard kinematic invariants denoted usually as 
t and u. We have added the squared screening momentum q'^ (see Sect. 
2) to the kinematic invariants in the denominators of Eqs. (^) and ( p!6| ) 
to account for the plasma screening effects. Furthermore, in Eq. (|l^) Civ 
and CiA are, respectively, the vector and the axial vector weak coupling 
constants associated with four-tail neutrino-emission vertices on electron 
(i = e) or proton (i = p) lines of the Feyman diagrams. For emission of Vei^e 
pairs at i = e (charged + neutral currents), one has CeV = 2sin^0w + 0.5 
and CeA = 0.5, while for emission of f^z^^ and i/^P^ (neutral currents only), 
Cgy = 2sin^0w — 0.5 and Cg^ = —0.5. In this case, is the Weinberg 
angle, sin^ ~ 0.23. Neutrino emission for i = p goes only via neutral 
currents and we have Cpv = (1 — 4sin^ ^w)/2 ~ 0.04 and CpA = 1.26/2 = 
0.63, independently of neutrino flavors. 

In our case, neutrino-pair momentum is sufficiently small {k ^ T <ti PFe), 
so that we can use the approximation of Festa & Ruderman [23] for the 
electron propagators: 
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and the nonrelativistic approximation for the proton propagators 

1 + 7° 



GAP2 + K) 



K) 



2uj 



(21) 



The standard tedious calculations using the Lenard identity [22] 



J dK J dkl {K-Ky- Kl) 



vr 



(22) 



yield 



Q 



12 



11 C12 y dpi ^ dp2 y dp'i ^ dp2 



12(27r)^ 

X Jl2/l/2(l-/0(l-/2) ■ 



For e-p scattering we have 
26C72 



2 

e+ "^p 



ep 



22 



2i? 



{PiK){P[K) 



(23) 



(24) 



where = Y^ui^lv + C^a)- Notice that Qep determined by Eqs. (H) 
and (|2^) appears to be very similar to the emissivity Qbr given by Eqs. 
(A7) and (A13) in the article by Haensel et al. [24] who considered neutrino 
bremsstrahlung due to Coulomb scattering of electrons by atomic nuclei 
in the NS crust. Equation (24) describes only neutrino emission associated 
with electron lines of the Feynman diagrams and does not include analogous 
emission associated with proton lines (the second term in Eq. (15)) which 
vanishes in the limit of nonrelativistic protons. 
In the case of ee scattering we have 



Cl+K\,h + J2-Jz) 



(25) 



where 
Jl = 
h = 



h , J2 
26 



(ill - 9s)(*22 



Y1 + Y2 + Z1 + Z2 , 
P,P2){P[P!,) + ^{P,P^){P2P[ 



X2 — I2 , J3 

3' D' 



+ 



(^^1^2 



^2 



{P2P[) 



{P[K){P!,K) {PiK){P2K) {PiK){P!,K) {P2K){P{K) 

Xi = 2' [{P,P2){PiP^) + iP2P[){PlP^)] 
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1 



it22-qW {PlK){P[K) 



+ 



Yi 



it 



11 



d)(ii2 
(PiP^) 



{P2P{] 



1 {P2P^) 

(tn-qW {P2K){PiK) 



{P2P^] 



{P{K){P^K) {P2K){P[K) {P2K){P^K) 
2\PiP2){P[P!,) 



{PiP[) 



{P2P[) 



{t22-qmi2- 

{PIP2) 

{PiK){P2K) {P^K){P[K) {P2K){P[K) 



(26) 



where Cg_ = ^^{Cly — C^^). The expressions for I2 and X2 are obtained 
from those for Ii and Xi by replacing 1^2. The expressions for Y2 and 
Z2 are obtained from those for Yi and Zi by replacing 1 jz^ 2 and 1' ^ 2' 
simultaneously. The quantities defined in Eqs. (|2^) possess numerous sym- 
metry properties. For instance, Xi is invariant with respect to replacing 
1^2 and 1' ^ 2', etc. The integrand in ( p3| ) can be shown to be highly 
symmetric as well. Using the symmetry properties one can easily prove that 
in Eq. (^) it is sufficient to use the expression 



In this case 



X 



Jee = 2 K\X -I^-Yi-Zi 



{PlP2){P{P^) + {P2P{){PlP^) {PlP[) 



it 



22 



{PiK){P[K) 



(27) 



(28) 



corresponds to "direct" ee scattering, i.e., to square of any single term 
(amplitude) in Ma or Mb in Eqs. (p^-(|l^); Ii corresponds to interference 
of amplitudes with exchanging electrons 1 and 2 (1 ^ 2, l' ^ 2'); while Yi 
and Zi correspond to interference of amplitudes with mutually transposed 
electron final states 1' ^ 2'. 

For further analysis, it is convenient to introduce an appropriate 4- 
momentum transfer Q = P2 — P2 = P\ — P[ — K = (Jl, q) from particle 
1 to particle 2 due to Coulomb interaction. Using the 4- momentum con- 



serving delta- function in Eq. (11) we can remove the integration over p'^ 
and over w, and then replace the integration over by the integration over 
q. It is convenient also to express Qep and Qee in the form suggested by 
Haensel et al. [24] for the emissivity Qbr of the neutrino bremsstrahlung due 
to scattering of electrons off atomic nuclei in the NS envelopes: 



Q 



12 



567 



(29) 
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In a Coulomb liquid of atomic nuclei, a dimensionless quantity L12 has 
meaning of a Coulomb logarithm. For neutrino ee and ep bremsstrahlung, 
we are interested in, L12 is a dimensionless function to be evaluated. General 
expressions for L12 are 

= /o ^wLg 3 /dpi /dp2 /dq /dk 



{q^ + q1Y {PlK){P[K) 



/i/2(l-/0(l-/2) , (30) 



Lee = ^ifi fnJ,fi / "^Pl / / dq / dk 



216^102^6^3^7 ^7 eie'ieae'a 

X i^' (X - /i - - Zi) /1/2 (1 - /0(1 - /2) • (31) 

To obtain Eq. (|30|) we set e2e2 = denominator of (^) because 

the protons are nonrelativistic. 

Since charged particles are strongly degenerate, the main contribution to 
L12 comes from those transitions, in which the particle momenta before and 
after collisions lie in the narrow thermal shells around their Fermi surfaces, 
|e — /i| < T. In Eqs. (^) and (^Tj) we may set dpj = m*pYi de^ dfij, where 
dilj is solid angle element in the direction of pj. We also put ei = pi in 
Eqs. ( ^0|) and (|3l|), and £2 = P2 in Eq- (^T|), respectively, since electrons are 
ultrarelativistic. 

Following Haensel et al. [24] we introduce the quantities q = q^ + q^ , 
k = ki + kr and p'/ = Pi — qt = Pi + q^ + k , where qt corresponds to purely 
elastic Coulomb scattering while q^ takes into account inelasticity; vector 
p'/ is directed along q^ but has the same length as pi, i.e., = |pi|; k^ 
and kj are the orthogonal vector components parallel and perpendicular to 
p'/, respectively. Notice that qt = 2pi sin(??/2), where i? is an angle between 
Pi and p'/. Strong electron degeneracy implies qr <^ PYe , k <^ p-pe and 
^ ^ PFei although qt can be comparable to ppe for large- angle electron 
scattering events. From geometrical consideration and energy-momentum 
conservation in analogy with expressions given by Haensel et al. [24] we 
obtain 

= (qt + ^.rf = qt +ql- qt (qr/pFe) , 

PiPi ~ {qt + + 2qt • k)/2 , 

PlK « p,{qr -n)-cifk, P[K « pi{qr - Q) , 

il, + U = El — e'l qr + kr , 

K^=u;^- kl - k^ ^ kl - k^ . (32) 

Here, k^ = {qr — Q){2u; — qr + and the condition > yields ko > kt 
and qr > u: > {qr — 0)/2 > 0. 
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In what follows, we consider ep and ee scattering separately. 



4. Neutrino ep bremsstrahlung 

Let us evaluate Q^p from Eqs. ( p9| ) and pO|): 

567m* r r r ^ r ^ f'^'^ ["^^^^ 

10 Jo 

dio I 

l{qr-n)/2 

{kl - kl) ql [1 - qVi^ple)] 



Ler> = j d.i j de, J dn, J dn, d^, dq, q, 

/OO roo rliT rko 

dqr / duj d^ dkt kt fl /2 (1 - /() (1 - /^) 

-oo J(qr-n)/2 Jo Jo 



eie[ {q^+q^,-n^ + qir (qr-rnqr-n-qt-k/pFe) ' ^ ' 

where we have introduced integration over dq = difg dqt qt dqr with cylin- 
drical axis along pi and integration over dk = dip kt dkt duj with cylindrical 
axis along p'(. In the latter integration we have used the equality dkr = div; 
ip is an azimuthal angle of k with respect to the p'^-axis. 

Consider the case of qs <C PFe typical for the NS cores, where qs is the 



plasma screening momentum (Sect. 2). Integration in Eq. (33) can be sim- 
plified because the main contribution comes from the values (/s < ^ PFe 
which correspond to the small-angle approximation. In this approximation, 
we can omit q^/{2pYef' in the numerator and qt • k/ppe in the denomina- 
tor of the integrand. Furthermore we may set l/{£ie'i) = l/ppg the 
integrand of ( |33|) and take into account that the energy transfer Q from an 
electron to a proton is 

O I ^ ( ^ P^^l Q (1A\ 

fi = e2-e2~— P2-q+— ~ — - cos 'dq2 , (34) 
m* \ 2 1 m* 

where t?g2 is an angle between q and p2. Then the integrand in Eq. ( p3[ ) 
depends only on the relative positions of p2 and q. 

First we can fix positions of pi, p2 and q and perform the following 
integrations: 



27r rko 

dip / dkt kt{kl - kl) = 27r(g,. - O)"' 
Jo 



UJ — -{qr — fl) 



(35) 



A sequence of angular integrations yields / dili / dipq J d^l2 ■ ■ ■ = 167r^ / 
dcos 'dq2 ■ ■ ■ Finally, we are left with the only angular integration in Eq. ( p3[ 



, that over dcos'!?g2- According to Eq. (|34[) , the variation range of cosT?q2 is 
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related to variation range of 0,. Let us consider two extreme cases. 

Low temperature limit (T <C qsPFe/i^^p ^ Qs) corresponds to thermal 
energies T much less than recoil energy of protons. This case is most im- 
portant for practical use [15]. Let us replace integration over cost?ij2 by 
integration over Q: cosT?q2 = ^p^/{QtPFe), where {qr — 2uj) <^<qr- The 
Pauli principle strongly restricts possible values of as well as final states 
of electrons and protons. 

The condition T <^ Qs leads to the inequalities qr <^ qt, uj ^ qt, ^ ^ qt- 
The integration over dqt in Eq. (|33|) yields 



2ppe fl2 1 / vr 

d'^tTI^^^ =^ TT— > (36) 



where fjs is given by Eq. (|8|). 

The integrations over ei and £2 are standard: 

/+00 r+oo 
dei / de2/i/2 (1-/0(1-/2) 
-00 J —00 

= J+uyrl 1 for l^>0,(. + f7)>0; 

= 1 -%\-l)/T for n<0,(^ + l^)<0, (37) 

where e'l = ei — lo — Q and £2 = £2 + ^- Now we are left with integrations 
over duj dfi dqr in six domains restricted by the inequalities qr > il. and 
1^ > {qr — 0,)/2 > (see Eqs. (^)). Three of them correspond to > 
[{a) n > 0; (b) n < 0, Lo + n > 0; (c) LJ + Q < 0], while three others 
correspond to > qr > [{d) u + Q > 0; (e) 2uj + > > u + (/) > 
2uj + ^} > qr]. Following Haensel et al. [24] we can introduce dimensionless 
variables u = qr/T, w = v = co/T, vq = {qr — 0,)/{2T). Six domains 

of integration over dv dw dvo in Eq. (^) are reduced to four ones, and 
integration over dvo is trivial (/q dvo {v — vq)'^ = 

Finally, using Eqs. (H), (|5|), (||), and (|3) we obtain 



l89m;'T' 3vr3 ( m;Ty 

Lep ~ „7 5 4 r] = — — ^ , 38 



where 



00 r roo 

dvv"^ < / dw 



w {v + w) 



(e^+«'-l)(l 
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/■^ , w (v + w) f" ^ w (v — w) 

+ / dw- — ■ V,,, ^ + dw- ^ ' 







(e^+«' - 1)(1 - e-"') Jo (e''-"' - l)(e"' - 1) 



oo 



+ / dw ^ , U 1646.7, (39) 

hv (1 -e''-"')(e"' - 1) J ' ^ ' 

as obtained by numerical integration, and ^ = 315 r//(167r^) ~ 3.417. 

Kaminker et al. [15] proposed to estimate Lep by rescaling the Coulomb 
logarithm Lbr calculated by Haensel et al. [24] for scattering of ultrarela- 
tivistic degenerate electrons by atomic nuclei in liquid state. In the context 
of ep scattering, Lbr is valid at Tpp ^ T <C Tpe. The rescaling procedure 

has been based on the similarity criterion: Lgp = -Lbf z^gp/ , where z^ep 
is an effective ep collision frequency for the conditions under consideration 
{T <^ Qs < Tpp <C Tpe) and vip^ is an effective frequency of elastic ep col- 
lisions at T ^ Qs (e.g., [25]). For T <C Qs, the authors used the collision 
frequency Vep that determines the electron thermal conductivity (e.g., [20]). 
One can easily verify that the result given by Eqs. ( |38| ) and ( |39[ ) differs 
from the crude estimate by Kaminker et al. [15] by a small factor given 
by (^). Thus Kaminker et al. [15] have strongly (typically, by two orders of 
magnitude) overestimated Qep, obtained in the present paper by an accurate 
evaluation. The nature of this overestimation comes from using the effective 
thermal-conduction collision frequency i^ep in the rescaling procedure. The 
thermal conduction at T <C (7s is associated with small momentum transfers 
{q ^ T). Corresponding collision frequency appears to be greatly enhanced 
as compared to the frequency of collisions with higher momentum transfers 
(q ^ Qs) which should have been used in rescaling. For instance, it would be 
appropriate to use the ep collision frequency which determines electron elec- 
tric conductivity (e.g., [26]). Therefore, the rescaling proposed in ref. [15] is 
basically correct if one uses physically justified effective collision frequencies. 
Finally, we obtain (in the standard physical units) 



_ TT^CGle^C!_^m;^ 

Qep — „ J.9 o 4 ^p {^bT) Rep 

945 Tl^c" Vs Ppe 
3.69 X 10^^ / m. 



2 / \ 1/3 

(40) 




Here we have introduced the factor Rep that describes suppression of the 
emissivity by the proton superfluidity. We introduced also Tg = T/10^ K. 

If protons are normal (T > T^p), one has Rep = 1 in Eq. (||). The 
factor Rep reduces neutrino emission at T < Tcp. It should be the same 
as the factor i?„p that describes reduction of the bremsstrahlung neutrino 
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emission in np collisions by proton superfluidity [3]: 
1 



2.732 



A = 0.9982 + J(0.0018)2 + (0.3815y)2 , 



Aexp (^1.306 - ^(1.306)2 + y2 
+ 1.732 5^ exp (3.303 - a/(3.303)2 + V 



B = 0.3949 + -y/ (0.6051)2 ^ (o.2666y)2 , (41) 

where y = A/T (Sect. 2). 

Moderate temperature limit {qsPFe/nT-p <C T <S Qs)- In this case thermal 
energy is much higher than the recoil energy of protons. Therefore, as 
follows from Eq. ([3^), T ^ Q. To estimate Lgp in this limit we can omit 
in comparison with all quantities ~ T in the integrand and in the integration 
limits in Eq. (^). Note that the condition kg ~ qr{2uj — Qr) > obviously 
gives Qr > and uj > qr/2. Then we easily perform integrations in Eq. ( p5[ ) 
and subsequent integrations over d^}2 dipq and dili. In the limit of $7 ^ 
standard integration yields 

ioT 



Jdei y'd£2/i/2(i-/0(i-/2) = ^jjrrrj- (42) 



which gives 



567m! r /-ZpFe qf roo roo v'^{v-u/2f 
' " ' au I av : — 




m/2 



e" 



(43) 



Here Lep is close to the Coulomb logarithm for neutrino bremsstrahlung due 
to collisions of electrons with nondegenerate protons in the low-temperature 
limit (Tpp < T Qs <^ ?Fe)- The factor 3m*T/p|p in Eq. ( |4^ ) appears due 
to degeneracy of protons in the case considered here. As a result, we obtain 
the emissivity (in the standard physical units) 



Qep = ion-fc9 8 2 V^bT) L Rf.p 

189n, c*ppg 



1.89 X 10^7 ( ^) (JJ^)'^' 4? ^epTj erg cm"^ s-\ (44) 
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Note that at low temperatures Qep has the temperature dependence Qep oc 
T^, which differs from the case of moderate temperatures, Qep oc T^, as well 
as from the case of nondegenerate protons, Qep oc (T > Tpp). 

The case <C T < Tpp has no importance for astrophysical applications. 
It can be shown that in this case the emissivity is again given by Eq. (E4|) 

with somewhat different (temperature dependent) Coulomb logarithm . 



5. Neutrino ee bremsstrahlung 



Let us consider now the neutrino emissivity Qee for ee bremsstrahlung. 
We restrict ourselves by the same condition Qs <C pfc (small-angle approxi- 
mation) as for ep bremsstrahlung. In contrast to the previous case, only the 
region of relatively small temperatures, T <C (7^, is important for practical 
applications. 



The integrand in Eq. (31) depends on the relative positions of pi, p2, 
q and k. In particular it depends on angle G between pi and p2- Thus we 
may direct pi along the z-axis and place p2 in the (x2;)-plane. Let (i?gi, ipq) 
be polar angles of q and {'&ki-,^k) be polar angles of k. Let us introduce 
also an angle i?g2 between q and p2, an angle t?fc2 between k and p2, and a 
relative azimuthal angle ipkq between k and q. 

We can use relations ( p2D as for ep bremsstrahlung. While evaluating 
the right hand side of Eq. (|3lD it is useful to introduce also an additional 
auxiliary integral / dr2(5(e2 — s'2 + Using the relation £5 = |P2 + q| ~ 
^2(1 + (q/pfc) cos'!?g2) we can present the 5-function under this integral in 
the form 

d{e2 -e'2 + n)^-6 (cos - -\ ■ (45) 



Equations (|31 



(PI), m, and (PD yield 



J^ee. — ^ee '^ee 



567 



dei y" de2 y"drJi y"d02 



ko 



r2pYE foo pQr /•27r pccj 

X / Mtqt I Mr dQ dip did I dkth 

JO JO J-00 JO J(qr-n)/2 JO 



X - I COSl/q2 — — 



(l-cosG)2 {ki-k^)q^ 



q \ ■ qJ {q! + q?-^' + ql? {qr - 

l + {Ce-/Ce+f (qr-n) 



1 



CO — k COS 1? 



k2 



^/i/2(l-/{)(l-/2) -(46) 
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Here we have introduced the same notations as in Eq. (^), and we have 

divided L^e into two parts, L^ee and L^ee ■ The first part, L^ee , corresponds 
to the first term (equal to 1) in square brackets which comes from X in 

Eq. (^) ("direct" ee scattering). The second part. Lie'*, corresponds to the 
second term which comes from Ii in Eq. (^) (interference of amphtudes 
with exchanging electrons 1 and 2). It can be shown that the terms Yi and 
Zi in Eq. (^) (interference of amplitudes with mutually transposed electron 
final states) contain small parameters ~ qr/pFe and ~ (qt/pFe)'^ and can be 
neglected. 

Let us use the relationship cos i9g2 = cos Q cos + sin G sin '&gi cos (pg in 
Eq. (^5|) and express cos ' dq2 in the (5-function through integration variables 
qt, qri ipq, and Q in Eq. (|4^). We also need to express cosT?fc2 through kr, 
kt, if, and 0. For this purpose we take into account additional approximate 
relations (justified at ppe ^ qt) 

cos 'i?g2 ~ {qr cos @ + qt sin 6 cos ^Pq)/q , 
cos ■i?yt2 ~ i^r COS Q + kt sin cos (pk)/k , 
cos ipk ~ sin ipkq ~ sin ip . (47) 

The third relationship is valid at qr <^ qs ^ qt, ^ ^ qt, and uj ^ qt. First we 

evaluate L^ee in Eq. ( ^). The integration over d^p is trivial, the integration 
over kt is reduced to (pSj), and the integration over dipq can be done with 
help of the delta-function 

d [cos'&q2 = — - d cos ipq — - — . (48) 

V q J qt smS \ qt sm 9 / 

Now we can integrate over df^i and dO,2 in a straightforward manner: 

[ [ dl7idl72(l -cose)2(sinG)"^ = 127r^ (49) 

The integrations over dei and de2 can be performed using Eq. (|37[) . We are 
left with a 4-dimensional integral (over do; dO dqr dqt). Using Eq. ( p6|) and 
introducing the same dimensionless variables as in Eqs. (^) and (]39|), we 
obtain 

,(11., 567 Sir' ( T \^ 



where r] is defined by Eq. ( p9| ) and = 189?7/(167r ) = 2.05. 

Making use of Eqs. (|4^) and performing straightforward but tedious cal- 



(2) 

culations we can also evaluate Lee • It appears to contain a small numerical 
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factor w 0.05 as compared to l'^J ■ Note that we have neglected other in- 



terference terms in Eq. (31) which could be of the same order of magnitude 



or even larger than Lee • Therefore we will use only the main term (|50|) for 



applications. 

Finally the neutrino emissivity Qee can be written in the form (in the 
standard units) 

0.69 X 10^^ fUeV^^ o 1 

« — r| erg cm"^ s~^ (51) 

Vs \no/ 

As in the case of Qep, this result differs from a crude estimate given by 
Kaminker et al. [15] by a small factor <C 1 due to an unfortunate use 
of inappropriate thermal-conduction ee collision frequency in the rescaling 
procedure. 

Let us emphasize that the emissivity Qep is reduced exponentially by 
strong proton superfluidity, while Qee is affected by the superfluidity in a 
much weaker manner, only through the plasma screening parameter (^). If 
protons are normal, they provide the major contribution into the plasma 
screening. If they are strongly superfluid, a weaker electron screening be- 
comes important, which enhances Qee, but not to a great extent (see below). 

Let us remind that the neutrino emissivity Qee for nondegenerate elec- 
trons in the neutron star crusts has been calculated by Cazzola & Sagion 
[16, 17]. We have compared their numerical results (Fig. 2 in ref. [17]) with 
our Eq. ( |5l| ) for several temperatures (Tg = 5, 7, 8, 10) at those values of 
P ^ 10^ g cm~^, where the relativistic electrons start to become degener- 
ate and both approaches are expected to be qualitatively valid. We have 
obtained reasonable agreement of the results. 

Finally let us make a few remarks about neutrino bremsstrahlung reac- 
tions involving muons, . Muons can appear in sufficiently dense matter 
of neutron star cores, provided the electron number density exceeds the 
threshold value Ue > {m^c)^ / {^i^'^fc') = 0.0052 fm~^, where is the muon 
mass. For realistic equations of state, this happens at p > (1.5-2) po- The 
muons, like electrons, form almost ideal degenerate gas, but they are mainly 
non-relativistic (near the threshold density) or mildly relativistic at much 
higher densities. Beta-equilibrium condition implies fi^ = fie, where fi^ is 
the muon chemical potential. The number density of muons is much smaller 
than the number density of neutrons or electrons for all densities. In muonic 
matter, several new neutrino bremsstrahlung reactions are allowed, result- 
ing from /ue, /u/x, fip scattering. The neutrino emissivities in these reactions 
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can be evaluated in the same manner as the emissivities due to ep and ee 
neutrino bremsstrahlung. 

For instance, in the nonrelativistic hmit (pp^t ^ ?7i^c) at T <C g^, in 
analogy with Eq. (Eof), we obtain 



4.68 X 10^2 fno\^/^ . n . 

' r| erg cm-3 s"^ , (52) 



Vs \n. 

where is the muon number density. Equation (|52| ) differs from Eq. ( ^0|) 
by a small factor (m^/m*)^. 

One can see that, under considered conditions, the neutrino emissivity 
Qeu is much smaller than the emissivity Qee due to ee scattering (cf. Eqs. 
(3 and (|l|)). 

Our estimates show that the neutrino bremsstrahlung reactions involving 
muons are less efficient than analogous reactions involving electrons, for all 
densities typical for neutron star cores. 



6. Discussion 



Let us discuss the efficiency of various neutrino emission mechanisms in 
superfluid NS cores. For illustration, we adopt a moderately stiff equation of 
state proposed by Prakash et al. [27] (the version with compression modulus 
Kq = 180 MeV and with the same simplified symmetry energy Sy which 
was used by Page & Applegate [28]). We assume that dense neutron-star 
matter consists of neutrons, protons and electrons (no muons and hyper- 
ons). The nucleon effective masses will be set equal to 0.7 of the masses of 
bare particles. The equation of state allows the direct Urea process to op- 
erate at p > 4.64/90. We consider two densities, p = 2po, as an example of 
the standard neutrino cooling (the direct Urea is forbidden), and p = 5po, 
representing the case of the cooling enhanced by the direct Urea process. 
Both nucleon species, neutrons and protons, are supposed to form superflu- 
ids. We consider a triplet-state neutron pairing ('^P2), with zero projection 
of the total momentum (mj = 0) of a Cooper pair onto quantization axis, 
and a singlet-state proton pairing (^Sq). Since the critical temperatures 
Ten and Tcp of both superfluids, as calculated from microscopic theories, 
are very model dependent (e.g., refs. [5, 6, 29-32]), we do not make any 
specific choice of microscopic superfluid model but treat Ten and Tcp as free 
parameters. 
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Figures 1 and 2 show temperature dependence of the neutrino emissivi- 
ties for various neutrino production mechanisms at p = 2pQ and p = 5po , 
respectively, assuming rather strong neutron and moderate proton super- 
fluids, Ten = 8 X 10^ K and T^p = 2.5 x 10^ K. 




IgT [K] 

Fig. 1. Temperature dependence of neutrino energy emission rates for various 
neutrino reactions in npe matter at p — 2po, Ten = 8 x 10^ K and T^p — 2.5 x 
10^ K. The neutrino reactions are: modified Urea (MU, sum of n and p branches), 
nucleon-nucleon bremsstrahlung (NN, sum of nn, np and pp branches), ee and ep 
bremsstrahlung. Cooper pairing of neutrons (CPn) and protons (CPp). 

In Fig. 1 we show neutrino emission from the modified Urea process (|l]) 
(sum of the proton and neutron branches) , nucleon-nucleon bremsstrahlung 
(|) (sum of nn, np and pp reactions), Cooper pairing of neutrons and protons 
(|), and ep and ee bremsstrahlung processes (|5|) and (^. The emissivities 
in the last two reactions are obtained in Sects. 4 and 5. The emissivities 
of the modified Urea and nucleon-nucleon bremsstrahlung reactions (||) and 
(^) are taken as from Levenfish & Yakovlev [7], with proper account for sup- 
pression of the reactions by the neutron and proton superfluidities [3] . The 
neutrino energy emission rate due to pairing of neutrons is plotted in ac- 
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cordance with the results by Yakovlev et al. [11, 12]. The neutrino emission 
due to pairing of protons is calculated taking into account the relativistic 
corrections [13] to the non-relativistic result [11, 12]. The allowance for 
the relativistic corrections increases production of neutrinos due to proton 
pairing by more than one order of magnitude. The neutrino emissivity due 
to ep bremsstrahlung is calculated in the low temperature limit (see Sect. 4) 
from Eqs. (^) and (^). We use Eq. ( [5l| ) to obtain the neutrino emission 
due to ee bremsstrahlung. 

At T > Ten the neutron-star matter is nonsuperfluid, so that the mod- 
ified Urea process and nucleon-nucleon bremsstrahlung play leading roles. It 
can be seen that the nucleon superfluidity reduces the Urea and bremsstrahlung 
neutrino reactions involving nucleons ((||), (|2|), and (^)) at T < 10^ K. On 
the other hand, the nucleon superfluidity produces two bumps of the par- 
tial neutrino emissivities at temperatures slightly below two critical tem- 
peratures. Ten and Tcp, due to Cooper pairing of neutrons and protons, 
respectively. The emission due to pairing of neutrons dominates over other 
mechanisms at T < 6 x 10^ K, while the emission due to pairing of pro- 
tons dominates at T < 1.5 x 10^ K. If the Cooper pairing emission were 
absent, the total neutrino emissivity at T < 10^ K would be 2-4 orders of 
magnitude smaller due to strong suppression by the nucleon superfluidity. 
The Cooper pairing can easily turn suppression into enhancement at not too 
low temperatures [11]. However, the Cooper-pairing neutrino reactions are 
eventually suppressed by superfluidity, when temperature falls down much 
below the critical temperature of neutrons or protons (T < 0.2Tc7v)- For 
T > 10^ K, the neutrino emission due to ep bremsstrahlung turns out to 
be more intense than that due to ee bremsstrahlung but both processes are 
much weaker than the other neutrino processes. With decreasing T, the situ- 
ation becomes drastically different. At T < 2.5 x 10^ K, ee bremsstrahlung 
becomes the leading mechanism of neutrino emission since it is the only 
mechanism which is almost insensitive to the superfluid state of matter and 
suffers no exponential suppression. Let us emphasize that its role is even 
more important for stronger superfluidities (for higher Ten and Tcp). 

In Fig. 2 we plot the contributions from the same neutrino reactions as in 
Fig. 1 and add the emissivity produced by the powerful direct Urea process 
(^). The reaction (^ is treated as described by Levenfish & Yakovlev [7], 
suppressed by nucleon superfluidity as in [21]. The direct Urea process is 
seen to dominate all other reactions at T > 2 x 10^ K. However, in spite of its 
high efficiency in nonsuperfluid matter, it is strongly suppressed at low T by 
nucleon superfluidity, so that ee bremsstrahlung becomes the leading neu- 
trino emission mechanism at T < 2x 10^ K. Therefore ee bremsstrahlung can 
play leading role even in the case of neutrino cooling in which the direct Urea 
process is open (by momentum conservation rules) but strongly reduced by 
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IgT [K] 

Fig. 2. Same as in Fig. 1 but for p = 5 po- Direct Urea (DU) process is operative, 
in addition to neutrino reactions shown in Fig. 1. 

superfluidity. Since the superfluidity reduces the direct Urea emissivity, it 
reduces also strong difference between the standard and enhanced cooUng. 
The temperature range of the ee bremsstrahlung domination is narrower for 
rapid cooling. 

Figures 3 and 4 display the domains of Ten and Tcp, where different neu- 
trino production mechanisms dominate at T = 2 x 10^ K for the standard 
and rapid cooling, respectively. Dashes show the lines T = T^n and T = T^p. 
In the left lower squares restricted by these lines, neutron-star matter is non- 
superfluid. Accordingly, the modified Urea process in Fig. 3 and the direct 
Urea process in Fig. 4 are the most powerful neutrino reactions in these do- 
mains. In the regions to the right of these domains neutrons are superfluid 
but protons are not, while in the regions above these domains protons are 
superfluid and neutrons are not. In the upper right squares enclosed by the 
dashed lines both nucleon species, n and p, are superfluid. One can see a 
variety of dominant mechanisms regulated by superfluidity. If both n and 
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7 8 9 10 

Ig Ten [K] 

Fig. 3. Domains of the neutron and proton critical temperatures Ten and Tcp in 
which different neutrino emission mechanisms dominate in a NS core at p = 2 po 
and T = 2 X 10^ K. Dashes show the hnes of T = Tc„ and T = T^p. 

p superfluidities are relatively strong (upper right corners of Figs. 3 and 
4), they switch off all the neutrino emission mechanisms involving nucleons, 
and therefore, the neutrino emission due to ee bremsstrahhmg dominates. 
This region becomes wider with decreasing temperature T (e.g., during NS 
cooling). Therefore, one of two reactions considered in this article, neutrino 
ee bremsstrahlung, can be the dominant neutrino reaction in both cases of 
standard and enhanced cooling. Wc will not discuss all variety of mecha- 
nisms which may dominate at different Ten, Tcp and T but notice that all of 
them are easily implemented into the cooling codes since their emissivities 
are described by simple fitting formulae [3, 7, 11, 12, 21]. The present article 
supplements this description with new neutrino reactions. 
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7 8 9 10 

Ig Ten [K] 

Fig. 4. Same as in Fig. 3 but for the neutrino emission enhanced by the direct 
Urea process at p = 5 po- 

7. Conclusions 

We have analyzed two neutrino emission mechanisms in the cores of 
neutron stars: neutrino ee and ep bremsstrahlung. We have obtained sim- 
ple expressions for calculating the neutrino energy emission rates in these 
reactions (Sects. 4 and 5). We have shown (Sect. 6) that the neutrino ee 
bremsstrahlung can dominate in superfluid cores of NSs in some domains 
of parameters of neutron-star matter. Therefore this reaction should be 
included in simulations of neutron star cooling. 

One may expect that ee bremsstrahlung will be the main neutrino 
production mechanism in highly superfluid cores (Tc„ ^ 3 x 10^ K and 
Tcp ^ 3 X 10^ K) of sufficiently old NSs (age t > 10^ yr) during the transi- 
tion from the neutrino cooling stage to the photon cooling stage. 
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